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ABSTRACT 

Ultraf iltration of Black liquor was carried out using an 
asymmetric membrane, in a stirred batch cel 1 , modified to work on a 
continuous mode. Polyethylene Glycol-6000 a standard macromolecule, 
was used to compare the results. The effects of operating variables 
viz. , differential pressure , bulk concentration and stirrer speed 
was studied in de^p^on permeate flux an observed rejection for 
both the solutes. 

5 

Osmotic pressure model was used to analyse the experimental 
results .Adsorptive resistance was taken into account 
quantitatively .The polarised layer resistance was found to be a 
function of membrane surface concent rat ion ,appl ied differential 
pressure and flow pattern. To relate these, the following 
correlation was considered and tested with the experimental 
results . 
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CHAPTER-1 

INTRODUCTION 

Ultraf iltration (UF) is a rate governed membrane separation 
process in which pressure is the driving force.lt has been gaining 
prominence as a practical process for concentration and 

purification of macromolecular and colloidal species in solu ion. 
Theoretically and now in practice UF thus offers an attractive 
alternative to a number of unit operations in food 

processing, chemical processing , pharmaceuti cal and medical 
industries for the concentration, purification and sterilisation of 
macromolecular and colloidal solutions on both laboratory and 
industrial scale. 

Ultrafiltration avoids many of the deleterious effects of 
conventional processes associated with phase change, high 
temperature and denaturing of components.lt also serves as 
a process for additional purification for concentrated 
products.GlimeniusC 1980) has presented the observat ions , vi ewpoints 
and the trends on UF for various industries. 

Pulp and paper industry is one of the major industries where 
UF finds extensive appl i cat i ons . It is essential to the economics 
of the kraft process of pulping that the spent chemicals from the 
process be recovered, regenerated and reused. The spent chemicals 
are seperated from the pulp at the washing stage along with wood 
residues, and is called Black Liquor (BL). Black liquor usually 
contains 98% of the total alkali charged to the digester. The 
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composition BL varies considerably in different paper 

mills, depending on the cooking liquor used and the method of 
cooking. host of the alkali in the BL is present as sodium 
carbonate or as organic sodium compounds very similar in chemical 
properties to sodium carbonate. The greater part of the organic 
matter removed from the wood in cooking is combined chemically 
with sodium hydroxide in the form of sodium salts of resinous and 
other organic acids .Appreciable amounts of organic sulphur 
compounds are present in association with sodium sulphite and the 
rest of alkali is present as free sodium hydroxide and sodium 
sulphide. There are other impurities too , such as 

lime, i r on-oxide, alumina, potash and sodium chloride. 

Black liquor recovery is an integral part of kraft 
process. The conventional BL regeneration process starts from 
concentrating the BL (10—15% total dissolved solids) by 
evaporating the water, then burning in a furnace to seperate 
inorganics and to recover heat. The inorganic ash as smelt is taken 
up in water , caustisized with lime, and reused in the cooking 
process. Make up sodium and sulphur are added to compensate for 
losses. The inorganics support combustion and provides steam for 
process use. 

The conventional process requires large capital investment 
and also destroys the valuable organics portion of BL. Further^ 
dilute BL directly obtained from pulp industries is of the order 
lOO cubic metre per ton of pulp (Gleminius, 1980) and hence it 
becomes necessary to recover water. Bigger paper mills have the 
potential to do so, but smaller mills simply discard the 
BL. Despite t$ its low concentration from smaller paper mills., 
disposal of BL is difficult because of high 
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BOD, colour , toxicity , etc above the standards. 

With rising energy costs , increasing capital investments and 
stricter environmental regulations , it has become more and more 
important to find efficient and alternative methods to concentrate 
and fractionate industrial effluents.UF is one such favourable 
process. 

Ultraf iltration finds major applications in pulp and paper 
industries particularly for the treatment of bleach plant 
effluerts and for spent sulphite liquor (G1 imenius , 1980) . 

Ultraf iltration of black liquor from pulp and paper mills has 
gained paramount importance not only as a waste treatment but also 
for recovery of valuable inorganic chemicals and this has become 
the controlling factor of the economy of the industry. The 
recovered chemicals are said to be valuable enough to give a short 
pay-back time for the UF plants . 

Gleminius( 1980) has reported the following advantages of UF 
to treat black liquor. 

( i ) the liquor is concentrated as much as 25 to 50 
times by UF. 

(ii) even a 50% rejection of solids results in atleast 
a 90% retention of colour. 

(iii) atleast 50% BOD and even more COD is retained. 

Civ) higher molecular weight 1 ignosul phates are 

retained. 

(v) less energy is consumed. 

All the above factors have initiated a resurgency in this 
area. However^ still little attention is being paid to the recovery 
of chemicals from kraf t BL using ultraf iltration. 



In earlier works (Mlahra and Bhattacharya, 1984 ;Rajnlah and 
Bhattacharya, 1985;Poddar et a 1.1969) attempts were made to analyse 
the electrodial is ive and ultraf i Itrative separation of BL. The 
present work has been undertaken to observe the ultraf i Itrative 
separation of BL and to study the limiting flux phenomena in 
UF.The work utilises BL as the observing system while a 

macromolecular solute (PEG-6000) being utilised as a standard 
solution for comparit ive analysis. The obj ect i ves of the study can 
be classified as : 

(i) to study the effect of independent variables 

( pressure, concentrat ion, t ime and flow velocity which are inherent 
to the system) on dependent variables (permeate flux and observed 
rejection) for UF of BL and compared to PEG-6000. 

(ii) to analyse the factors affecting the concentration 
polarisation and the resistance of the polarised layer to the 
f low, uti 1 ising an asymmetric membrane. 

( iii ) to correlate the mass transf er limiting phenomenon 
with the experimental variables. 



CHAPTER- 2 


LITERATURE REVIEU 

H.1 Development of membrane s 

Membrane separation processes have been in use on 
a small scale for over a century, but it is only within a period 
spanning the last three decades that research and development has 
advanced these processes to commercial signif icance.The pioneering 
membranes developed for membrane separation processes were the 
micro— filtration membranes cast in 1855. Later the first 

flat membranes were developed in 1872 . , however the 

development of membranes and studies on membrane properties did 
not catch up until and after the World War ll.Upto the sixties the 
membranes were classified as tortuos or sponge like until the 
advent of Nuclepore membranes from the General Electric 
Co r po rat i on , wh i ch revolutionised over the cellulosic 

membranes. These membranes were thinner, less porous and more 
flexible. The study of structure and shape of the membranes was 
made possible due to the invention of electron microscope. 

Early studies to reduce the porosity were highly 
empi ri cal .Researchers were studying the effects of concentration 
of membrane substance in casting, drying time required, solvent 
ratio, additives, etc. In 1957 it was found that cellulose acetate 
had the greatest permeability and 


rejection 



(Keating. , 1977 ) -The development of an excellent hypotheaia which 
could accomadate every class of membrane initiated the art of 
casting membranes. 

Contemporarily Loeb and Sourirajan (1958,1959,1960,1961) 
developed a new technique to cast Reverse Osmosis(RO) and Ultra 
Filtration(UF) membranes .These were able to reduce the resistance 
to the flow without affecting the other characteristics of the 
membranes. These membranes were integrally skinned anisotropic 
membranes and they initiated both commercial and academic 
interest, intitially in desalination and subsequently in various 
other membrane* processes. 

In an attempt to increase the active membrane surface area 
the hollow fiber membranes came into existence .DuPont and Dow 
Chemical developed a variety of such hollow fiber membranes -More 
recently (Rozelle &l <tl . , 197S) composite membranes heve been 
developed. These are structurally similar to Loeb— Sour i raj an 
membranes, but nevertheless they provide a free hand in tailoring 
various characteristics of membranes. 

Details about History and Development of membranes are 
discussed (Resting. , 1977; Schweitzer. , 197*1) . 

2.2 Ultrafiltration ; 

Ultratf iltrat ion (UF) through semi-permeable 
membranes has been recently evaluated as a most promising method 
for separation, purification, fractionation and concentration of 
substances like proteins and other macromolecules .This falls under 
the category of pressure driven membrane separation 
processes -There is considerable dispute as to the distinction 
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between Reverse Oenoais , Ultra Filtration and Micro Filtration 
(MF).The most useful definition of these processes is based on the 
smallest particle or molecule retained by the membrane .According 
to Porter( 1975, 1977) , the processes are defined by the rated pore 
size or the molecular— weight cut-off of the membranes as follows* 


O.OOOlpm 0.001pm 0.02pm 


10pm 


R O 


U F 


M F 


o 

1 A 


o 

10 A 


o 

200 A 


o 

100,000 A 


The characteristic features of this process (Dr iol i - , 1978) 
are followings 


< i ) 

<ii) 

(iii) 

(iv) 


It is athermal and can be carried out at 
constant temperature when necessary. 

No phase change is required. 

The operation needs low hydraulic pressure. 
The molecules are not denatured. 



(V) 

(vi) 



Chemical agents are not required generally. 

The process is economical. 

Ionic strength and p of the concentrate are 
maintained cosntant. 

In UF the dispersed phase C i .e. , solute) passes less readily 
through the membrane than the solvent. The reasons attributed for 
this may be classified as (Ferry ., 1936;Tuwiner 1962) : 

(i* The solute may be adsorbed on the surface of the 

filter and its por *s (primary adsorbtion). 

(ii> The solute may be retained within the pores 

(blocking) . 

(iii) It may be mechanically retained (seiving). 

One of the problems found in UF is a marked decline of the 
permeate flux with time. Albeit the rejection acheived could be 
quite high for even salts^the flux was too low to be of any pract- 
ical uti 1 i ty. .This is mostly attributed to the deposition of a 
concentrated layer of macro— solute molecules at the membane— 

solution interface;the layer acting as additional barrier to the 
flux of the solvent and micro solutes. This build up of retained 

solutes and particles on the membrane acts as a dynamic membrane 

which becomes the limiting resistance to the flow and offers a 
higher retention capability than the support itself .Blocking and 
primary adsorption may also synecgate . 

These phenomena have led to the coining of phrases 

"Concentration Polarisat ion"(Mi chaels . , 1968)-a concentrated 

boundary layer formed at the membrane surface, which reduces the 
flux-and "Foul ing"-the irreversible attachment of solute on the 
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membrane . Concentration polarisation nay be considered reversible 
including reversible deposit ion. This reversible deposition or 
reversible fouling may not be isolated ( Aimar , 1988) . A prolonged 
concentration polarised UF plant may lead to fouling phenomena. The 
primary factors for these two phenomena are the concentration of 
the f eed, operating differential pressure across the 
membrane, velocity of the feed and temperature .Various secondary 
factors like the residual charge on the membrane, ionic strength of 
the feed, etc. , are also found to affect. 

The literature available in concentration polarisation is 
vast and varied and dates back to 1966. Gill (1966) has done 
pioneering work starting with RO desalination. This was followed by 
Jonsson (1969) who tried to quantify concentrat ion polarisation 
from the governing laminar boundary natural convection 
solutions. Goldsmith (1971) used geometric properties of the 
membrane pores and macromolecules in the solution to quantify 
retention and mass transfer .Kozins ky and Lightfoot (1972) 
developed a model for a special case of UF through a rotating 
disk. 

Thomas (1973) concluded that in terms of a simple eddy diff- 
usion model boundary layer decreases with Schmidt number .Kozinsky 
&t al 1974) also attempted a similarity solution and found that 
concentration polarisation was essentially independent of suction 
velocity distribution. Porter (1972) and Shen ®£ oil. (1977) combined 
the effects of mass transfer by diffusion and formation of filter 
cake and concluded both to be interactive with each other .Bel luci 
et clI . ( 1975) and Swaminathan (1978) showed the results of 



protein UF considering gel polarisation model . Jonsson e-t otl.(1977) 
concluded the transition from laminar to turbulent region was in 
the range Reynolds No. 1000 to 2500 in a parallel plate channel. An 
integral solution for the limiting case fo linear osmotic pressure 
was obtained (Probstein, 1979) .Jonsson (1980) has outlined the app- 
licabilty of various models for the transport of water and solute 
in UF/RO membranes. Guantitave models for membrane fouling have also 
been reported by Green *t «ai.(1980). 

The results of so’ ute mass balance solved by an integral 
method were shown to contradict the film theory model (Trettin 
a.1 ., 1980? 1981 ) (the film theory model supports that the permeate 
flux varies logarithmically with solute concentration.). 

Modelling for colloidal fouling was developed by Schippers 
clI .(1981). This was based upon the cake or gel filtration 
theory.lt was suggested that at the initial period of UF blocking 
filtration may take place insead of gel filtration. 

While in these experiments the anisotropic membrane was 
considered to have a uniformly distrubuted permeabi 1 i ty , Fane et 
<al.(1981) took regions of differing permeability into 
consideration for analysis. The gel polarised UF membrane was 
viewed as a physical analogue of a composite RO membrane .This gel 
polarised UF was said to depend on membrane permeability and 
surface properties.A correction factor was introduced for the 
available free area by them. 

The effect of secondary factors like solute size, solute hyd- 
rophobicity, p H and the ionic strength of the feed, etc., was studi- 
ed by Zeman ( 1983) .Multiple solute UF was found to show different 



results compared to single solute UF on the time of adsorption . 

The influence of solute properties and membrane surface 
characteristics on adsorption kinetics and hydronamics shows 
atleast two distinct controlling steps of adsorption 
(Mathiasson, 1983) .A strong variation of t ransmembrane pressure was 
modelled to quantify the degree of fouling (Flemmer, 1983) .Buckley 
(1983) developed a model to predict flux decline before and after 
the onset of gel polarisation. 

A summary of the concept of concentration 
polarisation, 1 imitations and the ways to reduce concent rat ion 
polar i sat iom, membrane fouling, solute membrane interactions 
etc. have been discussed in detail (Belfort <ei <atl.1983). 

The dynamics of polarisation in stirred and unstirred UF was 
studied by Chudacek et al .( 1984) .The simple Carman— Kozeny equation 
was shown to be a good aprroximatidn to the amount of solute in 
the polarised layer. Filtration theory was utilised to evaluate the 
specific resistance .Marked affection of water permeation and 
permeability due to adsorption of hydrocarbon impurities and trace 
inorganic ions also affects concentration polarisation and fouling 
(Errede, 1984) . 

Mi j mans et al. (1984,1985) have explained the flux 
limitations in UF in terms of osmotic pressure. On the other hand 
Jonsson (1984) suggests that a degree of osmotic limitation has to 
be considered along with gel polarisation model .Dickson at 
<al.(1984) and Do at ctl . ( 1984 ) have also discussed in this line. The 
limiting flux is independent of applied pressure and is strongly 



dependent on stirred and bulk concentration. 

Shu-Sen Wang «•(. a. I .(1985) have used least squares approxim- 
ation to determine the rate constant of gel growth. By defining the 
concentration profiles in the boundary layer they have modelled to 
determine the period of cleaning of the membrane surface.By utili- 
sing permeability data from sedimentation experiments, Wijmans at 
otl.(1985) have showed that osmotic pressure model and the boundary 
layer resistance model are equivalent. 

There is & particular interest in UF of Bovine Serum Albumin 

(BSA) and Poly Acrylic Acid (PAA). On these systems ,Choe 

oil .(1986) have conducted batch UF experiments. It was concluded 

that the polymer chain flexibility affects the fouling to a great 
H 

extent. p and ionic concentration affects the fouling of PAA but 
does not affect fouling of BSA. They also suggest that while for 
PAA the fouling is mainly due to pore plugging, BSA fouling is due 
to solute— membrane interact ion .The same result was obtained by 
Aimar ctl . ( 1986) .Adsorption of BSA was reported to follow 

Freundl i ch law. It was also reported that drastic changes in the 
filtration property of BSA due to solute-membrane interaction. By 
assuming the polarised layer to be less porous and highly resisti- 
ve to the f low,Suki <21.(1986) have derived a model. This model 
is derived from flocculation theory and it shows a high solute- 
membrane interaction. 

Nakao «t <21.(1981,1986) considered cake filtration theory to 
analyse concentration polarisation and fouling in unstirred 
conditions .This was followed soon by Aimar et ctl . (1986,1988) .In 
the later work they have put forward four models for evaluating 
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concentration polarisation. The convective deposition with constant 
removal of flux developed by them matches closely with their 
experimental results. A model based on adsorption isotherm and 
another on gel— layer growth have also been derived and discussed. 

There are a few literature review papers available for the 
mass transfer phenomena in UF.Of these Gekas *£ ctl.(1987) presents 
a detailed discussion on factors like porosity, wall 

roughness, change in viscosity and diffusion coefficient .Gill 
( 1988) has discussed the effect of viscosity on fouling and mass 
transfer using an integral model. The optimal operating conditions 
to be chosen are discussed in detai 1 , although there is a defect 
in the model due the assumption of constancy in permeate flux. More 
extensive review s are also available (Porter : * . , 1986;Phillip 

clI. , 1986;Muhammed, 1982;Belfort , 1984;Mazid, 1988) . 



CHAPTER-3 


THEORETICAL CONSIDERATIONS 


3.1 Flux through the membrane : 

In the absence of a solute , the flux , through a 
porous membrane is defined by Darcy's law, which states that the 
pure solvent flux is directly proportional to the applied pressure 
differential ,AP , and inversely proportional to the pure solvent 
viscosity, (Nguyen et al.,1979)s 


J = K AP _ ; K 
w appl 


R ^ 
m w 


C3.1) 


where R m ,is the membrane intrinsic hydraulic resistance.lt is a 
function of pore size, tortuosity, membrane thickness and porosity 
as follows: 


Sr <5 


R 


& a 


(3.2) 


Since UF membranes are plastic and can yield ( compact Junder 

pressure, R is also a function of the pressure history. Initial 
m 

membrane compaction occurs rapidly during start-up, whereas long 
term compaction occurs slowly aver the operating life of the 
membrane (K1 inkowski , 1984) .However , improvements made with thin-film 
composite membranes lessen this effect. 

Equation 3.1 no longer holds true when the solution contains 
retained solutes. This flux (J uf ) in an ultrafiltration run is 
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generally lover than the pure solvent flux In absence of solute 

» J w »and this is also decreasing with time. For example Fig. 3.1 

shows the variations of f lux ,J , versus the applied pressure 

differential ,AP ,for an UF run. 

appl 

In Fig 3.1 curve *a' represents Darcy's law for the pure 
solvent (J w ), while curve *b' represents the actual experimental UF 
flux (J ^J.From the figure it can also be noted that the slope of 
J u ^ is less than that of Darcy's law, whi le at higher pressures the 
flux, J u |fF becomes asymptotic to X-axis , thereby implying to be 

independent of applied pressure.lt may be further said that point 
’ c* is time dependent while point ' d* is time independent .These 

points * c'and *d' may be drawn as a function of time as on Fig. 3. 2 

If after an UF run the membrane is thoroughly washed with 

f 

appropriate washing solution and the solvent flux (J tf ) is 

determined under the same applied pressurse Ap appl » ® a y be found 

to be less than J w ,but still greater than J uf .The difference 

» 

between J and J is accounted for the irreversible fouling of the 
w w 

membrane due to the attachment of matter on the surface and inside 
the pores of membranes .This loss of flux can be visualised as 
additional resistance to the flux,R a . Hence 



(3.3) 


The difference between J uf and J w arises from a phenomena^ 
which disappear when the UF pressure is released or when the 
membrane is rinsed. There are two different types of phenomena that 




n 



Fig. 3*2 A graphical illustration of the possible ways for 
the drop in permeate f lux 
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may explain this: 


<i) the osmotic pressure model says a negative and 

reversible osmotic pressure is developed (Clifton 
et al. , 1984; Vi Ike r et ctl. ,1981), or 
(ii) as per the gel-polarisation model a depositive 

layer resistive to the flow is said to be formed 
which is reversible in nature. 

It has been impossible until now to distinguish these two 

phenomena, which get modified as soon as the UF mass transfer is 

changed in any manner .As a consequence all these phenomena may be 

characterized as if they reduce the applied pressure AP -,,by an 

appi 

equivalent osmotic pressure , ATI, down to an effective 

pressure ,AP^^ , defined as ( Aimar , 1988) : 

AP eff " AP ap P r An < 3 - 4 ' 


In recent times Nakao et ctl. (1985) have showed that both the 
above models are equivalent to each other, while Jonsson (1984) 
suggests that a degree of osmotic limitation along with 
gel-polarisation to be considered.The above two models may be 
theor it i cal ly represented as follows: 


J 


uf 


ap _ - An 

appl 

fj Tr +r +r 1 

si m p aj 


(3.5) 


where is the viscosity of the feed solution and ATI is the 
osmotic pressure difference between the feed and the permeate 
membrane surface concent rat ions . 

Experimentally R and R could be obtained by measuring pure 

** HI 
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solvent flux J end J from Eqns .3.1 An 3.3. 
w w 

The time dependent case of Eqn 3.5 can be modified as: 


J uf <l) 


ap . - An 
appl 


( R 


+R ( t )+R 


SL * P 


al 


(3.6) 


From Fig. 3. 2 it should also be noted that reaches an 
almost constant final flux (F) .The time coresponding to this 
J uf (F) is T(F) .This factor implies that at . J uf (F) ,R p ( t ) 
becomes constant ,R p (F> .Hence, Eqn. 3. 6 may be modified at this final 
level as: 


J 


uf 


(F) 


AP 


appl 


- An 


v Tr +r CF)+R 1 
H ® P aj 


(3.7) 


3.2 Concentration polarisation : 

In membrane seperation processes a phenomenon 
called concentration polarisation always exists , due to which the 
concentration at the membrane solution interface is always higher 
than the bulk solution concentration. This phenomena may be 
mathematically modelled by writing the solute mass balance for an 
infinitesimal element in the boundary layer and further 
integrating the following with the appropriate boundary 
conditions : 


J C - J C 

P 



(3.8) 




On integration v* obtain: 



J 


k In 


( 



(3.9) 


D 

where k = is the mass transfer coefficient. 

& 

For a stirred batch cell the mass transfer coefficient can 
be found by the following empirical correlation 

(Sherwood, T.K. , 1975) . 

^ 0.33 ^ Z 0.8 

X = 0.0443 (— ) {^) ( — <3 - 10 ' 

where , r = radius of the membrane (cm), 

v — kinematic viscosity of the feed solution (St), 

= angular velocity (rad/s) and 

2 

D = dif f usi vity (cm /s). 

Further the dif f usivity has been re lated empirically to 
molecular weight of polymeric solutions (Sherwood " .,1975) as 

follows: 

D = 2.74 x 10 _5 ( M ) _1/3 (3.11) 

3.3 Osmotic pressure model and Flory * s eq uation : 

Because of the concentration polarisation phenomena, the 
concentration at the membrane solution interface is always higher 
than that at the bulk. Due to this a negative osmotic pressure is 
developed which acts against the applied pressure, Hence the 
effective pressure differential is less than - * the applied 

pressure (Eqn.3.4) .In osmotic pressure model, this effective 
pressure dif ferential is used to calculate Jas follows: 



*■» Jt® 


J = 


AP 


appl 


An 


ij r r +r i 

8 1 n a J 


(3.12) 


Here All is the osmotic pressure difference between the 
membrane solution interface at the concentrate and permeate sides 
respectively. 


The osmotic pressure depends on the 

concentration 

and the 

molecular wei ht of 

the 

solute. From thermodynamics we 

can 

say at 

osmotic equilibrium 

the 

chemical potential 

of the 

solvent must 

have the same value 

- 

I on 

each 

of the 

membrane .Hence 






n 


(Hi - H°> 



(3.13) 
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v. R T 





where ^ is the chemical potential of the solvent in solution,/^ 
is the chemical potential of pqjje solvent at the same temperature. 

Me can express the chemical potential jMj as 
follows (Flory. , 1953) . 

F*i - = R J ( i" + ci - -j- i r z + x 1 y\ ] <3 - 14) 

From Eqns. 3.13 and 3.14 


n = 


RT 


V. 


( ln r i 


+ ci - 


□ r . 


x y 2 

*1 *2 


(3.15) 



where 


X 


volume of polymer molecule 


volume of solvent molecule 

v M / N 
p ' av 

V 1 / N av 
V P M 

(3.16) 


Here V, is the molar volume of the solvent , v is the 
1 P 

specific volume of polymer , and M is the molecular weight of the 
solvent .Eqn .3. 16 can also be written as 


X 


v M 
P 


v,M, 
1 1 


(3.17) 


where v^ is the specific volume of solvent and is the 

molecular weight of solvent. 

Y ^ and j'g of Eqn .3. 14 depend strongly on concent rat ion and 


can be expressed as: 


r 2 = G v p = C / p p 


(3.18) 


where C is the concentration, p is the density of the polymer at 

P 

the amorphous state ( = 1.125 gm/cc) (Vankrevelen. , 1972) , and 


Y ^ - 1 ~ (3.19) 

X . = 0.45 for PEG-water system (Vankrevelen. , 1972) . 

1 

So one can calculate the osmotic pressure from Eqns.3.14 


through 3.19. 
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3.4 Obaever fed and Real Rejection : 

Rejection is a measure of the ability of the membrane 

to retain a particular solute. Two types of rejection are generally 

used to describe the characteristics of a membrane , name ly the 

observed and the real rejections. 

The observed rejection is expressed as 

R = 1 - C / (3.20) 

o p • b 

wh~re 0^ is the permeate concentration and is '.he bulk 

concentration.But due to the concentration polarisation the 

membrane solution interface concentration is higher than that at 
the bulk. So the real rejection is 


R = 1 - C / C (3.21) 

r p ' m 

where C is the membrane surface concentration. 

11 

3.5 Velocity variation techniq ues 

Mass transfer coefficient (k) can be determined by 
using velocity variation technique (Nakao e-i ctl . , 1981 ) ,f or a 
solute of unknown physical properties .From Eqn.3. 10, it is known at 
any constant concentration, k is proportional to the stirrer speed 
(<o). So we can write, 

k = e (<o)®*® ; e is a constant (3.22) 

Using Eqns.3.20 and 3. 21, Eqn.3. 9 can be expressed as 


( 



= In 



J 
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] 


In 


(3.23) 



AS 


Subsltuting *qn. 3.22 into Eqn.3.23, 

1 - R q 1 - R r J 

ln ( _ “ ) = *" ( ] + V ' ( ~^T ) 

o r 


Equation 3.24 is a straight line plot of ln 



(3.24) 


versus (j/« 0 * 8 ) .The value of e and therefore mass transfer 
coefficient can be found from the slope of this plot. 
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EXPERIMENTAL WORK 


4.1 Instruments and Materials* 

4.1.1 Membrane s 
Type 

Nature 

Diameter 

M.W.C.O. 

Composition 

|J 

p range 

Maximum temperature 

Support 

Skin 

Sterilisation 

Suggested cleaning 

Storage 


Spectra/Py r » C5 

SPECTRUM* Medical Industries 
Inc»U.S.A. 

Flat disc, asymmetric, 
anisotropic and hydrophillic. 

76 mm. 

5000, sharp. 

cellulose acetate material 
or hydrol ised cellulose 

acetate . 

2 to 10. 

90' C. 

highly porous, 150 to 300 pm. 
0.1 to 0.5 pm. 

by formalin or hydrogen 
peroxide or ethylene oxide, 
by detergent with enzymes or 
2.57. hydrogen peroxide. 

1 to 2 % formaldehyde. 



4 . 1.2 


Stirrer : 



Diameter of the 
stirrer 

Height of the stirre 
above the membrane 

4.1.3 Stirrer motor s 
Ty » e 
Supply 
Horse power 
rpm 

Amperage 

4.1.4 Pump * 

Type 
Supply 
Amperage 
rpm 

Wattage 

4.1.6 Stroboscope s 

Type 

Range 

rpm per division 


56 mm. 

30 mm. 

shunt , continuous rating- 
220 V, DC. 

0.125 HP. 

1450. 

0.85 A. 

reciprocating. 

220V, 1Ph,50Hz. 

2.6A. 

1425. 

180W. 

1531— A 

General Radio Company, 
U.S.A. 

110 to 690. 

5 



4.1.7a 


Refract oneter: 



Model Pzo,Warsawa,RL1 .Poland. 

Range of measuring 

scale 1.3 to 1.7 

Accuracy 0.0002 1 .30<R.I<1 .42 

0.0004 1 .42<R.I<1 .70 

4.1.7b Weighing balance : 

Model 

Accuracy 

Maximum weight 

4.2 Experimental device s 

4.2.1 Apparatus : 

The main aim of the experiments were to examine 
the factors affecting the decrease in flux with respect to time. So 
it was nessary to conduct the experiments on a continuous mode and 
at steady state. Since the experiments were conducted using a batch 
cell ultrafilter the permeate had to be pumped back into the cell 
to maintain the bulk concentration constant .This would enable to 
conduct the experiment at steady state. 

In practice to recycle the permeate and there by maintain 
the bulk concentration constant was achieved using a high pressure 
metering pump. The metering pump itself was used to pressurise the 
system, thereby eliminating the need of a compressor or compressed 
nitrogen. The experimental set-up is shown in Fig. 4.1. 

The specifications of the cell (Fig .4.2) , designed and 


Sarat o r i us , 2004MP . 

O.OOOOIg 

166. Og 


fabricated were as follows: 



Stirrer 

motor 



F i g.4-1 Schematic diagram tor the experimental setup 
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1. Membrane, 2. Porous Polyethylene Support, 3/0'- Ring 
4. Cell Body , 5. Pressure Line Connector, 6. Pressure Gouge 
7. High Pressure Proved Stuffein Box, 8. Adaptor, 9. Stirrer 
10. Sample Outlet, 11. Feed Sample Inlet. 


Fig. 4.2 Ullraf i It rat ion cell. 


Material of 


31 

construction SS-316 

Useful volume 450 cc 

Residual volume nil 

Membrane diameter 76mm 

Effective filtration 

2 

diameter 26.4 cm 

Maximum testing 

pressure 500 psig 

Stirring faci 1 ity mechanical 

4.2.2 Solutions ^ 

Solutions of Poly Ethylene Glycol (PEG-6000) of 
molecular weight range 6000-7500 and Black Liquor(BL) were 
prepared using distilled water .PEG-6000 of A. R. grade was obtained 
from Fluka, England, whi le BL (approximately 16% total dissolved 
solids) was obtained from Central pulp mills, Surat. Solutions of 1% 
,2% ,5% concentrations of PEG— 6000 and BL were prepared. 

4.2.3 Design of experiments • 

Experiments were designed so that the effects of 
three major independent variables , pressure differential , bulk 
concentration and stirrer speed were studied. The experiments were 
done so that when two variables (e.g. pressure and concentration) 
were fixed, the third variable (e.g. stirrer speed) was changed for 
each experiment until all its levels were exhausted. 

The levels for concentration were 1% ,2% and 5% while the 
pressures were varied as 40psi ,80psi , lOOps i and stirrer speeds 
were varied as 325rpm,450rpm and 530 r pm. The standard levels of the 
variables were chosen as 80psi ,450rpm and 2% concentration .Of 
these^two were variables fixed while the third was different for 



each experiment. 



4.3 Experimental procedure s 

After the experimental levels of the variables were 
fixed, the experiments were conducted for each set of variables.The 
experimental procedure undertaken was as follows: 

The fresh porous membrane was placed on a porous support at 
the bottom of the cell. The cell was then assembled. Fi rst and 
foremost was to allow the compaction of the membrane by 
pressurising with distilled water at 120psi for 2 hours. The 
constancy of water flux beyond this time ensures no further 
compaction of the membrane .This reading is recorded for the 
membrane's solvent permeability. Once the cell was assembled with 
a well cleaned and rinsed membrane, the stirrer speed (for the 
required rpm) was adjusted by regulating the voltage supply 
through a variate. The rpm was measured by a stroboscope. The feed 
solution was pumped into the cell by the metering pump. Once the 
required pressure was attained the pumping of the feed solution 
was stopped. The permeate was recycled. The constancy of the 
pressure was adjusted by regulating the pumping rate of the 
metering pump. The permeate solution that was coming out through 
the outlet was collected and the time required for collection of 
lOcc of permeate was also noted. Flux was then calculated from this 
reading. Such readings were taken after every 15 minutes of the run 
until flux value was constant for three consecutive readings. All 
these collected permeate solutions were either recycled or 
replaced with freshly prepared solutions of equal 
concentrations, to the cell. There was a marginal increase of the 


feed temperature due to pumping of the solution, however for 
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practical purposes it waa neglected and aaauaed to be conatant 
(ambient temperature). 

After each run the whole cell was rinsed with distilled 
water thoroughly .The membrane was thoroughly washed and rinsed 
with distilled water for two hours to remove any deposit ion. The 
water flux was again checked at standard conditions and 

recorded.This procedure was repeated for each set of experimental 
conditions . 

4.3.1 Measurement of concentration s 

Refractive index calibration method was used to 
determine the concentration of PEG. 

To find the concentration of black liquor solution, lOcc 
of the solution was evaporated in an oven at 110’C in a pre- 
weighed pe tri-dish. The mass of the remaining solids after 
evaporation was used to calculate the concentration (total 
dissolved solids). 


4.3.2 Measurement of Viscosity and Density of Solutions s 

Viscosities were determined by an Ostwald Viscometer for 
various PEG and BL concentrations in a constant temperature bath. 

A correlation for viscosity and concentration was developed 
from these data (Appendix A-1). 

The densities were measured by using a specific gravity 


bottle (Appendix A-2). 



CHAPTER-5 


RESULTS AND DISCUSSION 
5.1 Selection of Black Liquor Concentration : 

The concentration and composition of black liquor (BL> 
varies , depending upon type of alkaline pulping, capacity of the 
industry and type of raw materials used for pulp 
product ion , etc .The bigger paper mills produces aruond 15% BL 
concentration (t.d.s) from brown stock washers , whereas small scale 
pulp industries discharge BL at a range from 1 to S% 
cocent rat ion.The recovery of BL from bigger paper mills is carried 
out with conventional technique (Chapter 1) whereas smaller paper 
mills simply discard. There is an urgency to find an economical and 
simplified procedure to recover the valuable organic as well as 
inorganic products from BL with simultaneous disposal of treated 
effluent to river st reams , under specified environmental rules. As 
highlighted in Chapter 1 ,UF is one such favourable technology to 
treat BL. 

The selection of BL concentrat ion for exper imentat ion plays 
an important role while studying flux and rejection behaviour of 
an asymmetric membrane. The following factors have to be considered 
before choosing the concent ration range for the experiments s 

(i) The activity and hence the osmotic pressure of the BL 
solution increases with an increase in the concentration of BL 
(Barrow, 1966) .This results in a sharp decrease in flux. The flux 
through the membrane becomes too low for any compar itive study at 
higher concentrations. 



(ii) Considering the batch cell with the stirrer designed 
for experimentation, it was observed proper to handle BL of 
concentration upto 5% only, as BL viscosity increases beyond this 
limit of concentration. 

(iii) Black liquor is basically composed of macromolecules of 
various sizes in alkaline solution. At slightly higher 
concentration, clogging of pores of the membrane takes place and 
reuse of the membranes becomes more difficult. 

For tht >e reasons, use of dilute solutions of BL is always 
recommended for experimental studies .Hence for this work 
concentrations of 0.5% to 5% BL were chosen to study the effect of 
operating variables on limiting flux and rejection. 

5.2 Selection of PEG Concentration : 

In order to make a comparit ive study with BL an unknown 
molecular weight it is essential that the concentration of 
standard solute is also of the same range(0.5 to 5%) . 

5.3 Effect of Operating Variables on Permeate Flux : 

Differential pressure , bulk concentration and stirrer speed - 
were the three variables varied to study the effect of flux as a 
function of time. 

5-3.1 Effect of Pressure on Permeate Flux : 

Fig .5. 1 shows the variation of permeate flux with respect to 
time for a 2% solution. The nature of the curves coincides with the 
trend shown in Fig. 3. 2, as expected. The flux for both PEG and BL d- 
rops initially and later attains an almost constant level. The flux 
at 2. 8x10”* and 6. < ?x1G^N/m*~ for PEG is less than that of correspo- 
nding BL flux. This suggests that although the bulk solution conce- 
ntration of both are same (2%), the average molecular weight of, t- 
he BL may be less compared to PEG. Considering this explanation of 



(s 01/01 oix) t n r'xn]j 




average molecular weight of BL to be leas than that of PEG, the r- 
ejection of BL solutes should be much less compared to PEG. It was 
indeed observed to be so and would be shown later . 

On the other hand for the intermediate pressure 
5 2 

(6.5x10 N/m ) the flux of PEG is marginally higher at the later 
stages inspite of being too low initially than BL.This may be 
attributed to the fact that BL being a polydispersed solute of 
wide range of molecular weight , some pores might have been blocked 
par r ial ly. Hence , after some time, one observes lower flux t tan PEG. 

The other fact to note is the time taken to reach the 
limiting flux increases with increase in pressure f o BL.Contrari ly 
this time to reach 1 imit ing flux for PEG is nearly same for all 
the three pressures. This implies that an increase in pressure not 
only increases the flux but also takes longer time to reach the 
constant flux for BL.But for PEG it only acts - to 
increase the f lux. Fur the r , with the increase in pressure there is a 
late development of concentrated boundary layer in the case of BL, 
whereas ,PEG solution does not show development of appreciable 
boundary layer as a function of time. This may be due to the low 
concentration of PEG. Black liquor, a peculiar polydispersed 
solution, builds up layer even at low concentration, thus showing 
different trend at different pressures. 

5-3.2 Effect of Stirrer Speed on Permeate Flux: 

It was obvious from Fig. 5. 2 that the permeate flux increases 
with increase in st irrer speed.lt can also be noted that the flux 
of BL is higher than that of PEG. On the other hand though the 
initial BL fluxes are higher for 47.12 rad/ s and 34.0 rad/s, the 
limiting flux of BL is less than that of the corresponding PEG 
curve .This was already discussed in Sec. 5. 3. 1, due 


to 


the 



Flux.J.,4 (xlCf m/m s) 



Fig. 5*2 Effect of stirrer speed on permeate flux 
for 2 °/o solution at A Pappl z 5*5 x 10 N/m 



polydispersed nature of BL. 

The time for reaching the limiting flux value for BL is 
distinctly decreasing with an increase in stirrer speed. This may 
be due to the possibility that the rate of deposition of polarised 
layer on membrane surface is faster, as Reynolds No. 
increases. Higher Reynolds No. increases turbulence within the cell 
which restricts the development of concentrated boundary 
layer, giving lower diffusive flux than convective flux 
(Fig. 3. 3) -Higher convective f lux , increases ;he possiblity of 
faster deposition of polarised layer on the membrane 
surface , because of the polydispersed nature of BL having tendency 
to precipitate on the surface of the membrane .This gives rise to 
an earlier break in the curve of Fig. 5. 2. 

It was observed that for PEG solution, this kind of trend is 
not significant. This seems to be logical considering the nature of 
PEG solute which does not precipitate. Hence at higher turbulence 
with more convective flux than diffusive and having much less 
deposition of polarised layer,the only possibility is to have 
steady permeate flux. 

5.3.3 Effect of Concentration on Permeate Flux : 

It is a well known fact that as the bulk concentration of 
the solution increases the permeate flux decreases .This is shown 
in Fig. 5. 3. Higher concentrations realises higher osmotic 
pr e s s ur e . Hen ce lower flux is observed. Further , higher 
concentrations may be giving higher polarised layer resistence 
which is inversely proportional to f lux. However , the trend of 
curves are identical as bulk concent rat ions increase at a fixed 
pressure and stirrer speed. As for PEG solutions, it is apparent 
(Fig. 5. 3) that there is not any significant polarisation as a 
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function of tla*,aa we increase the bulk concentrations. 


5-4 Effect of Operating Variables on Observed Rejection s 

The ef feet of operating variables on observed rejection is 
shown in Fig. 5. 4. The permeate concentrations did not vary in a 
single run for either PEG or BL. Observed rejections were 
calculated for each run from Eqn.3.20. 

5.4.1 Effect of Pressure on Observed Rejection s 

The curves obtained for PEG ^nd BL were in contrast. For PEG 
the observed rejection was independent of pressure at low 
pressures ,whi le at high pressures it decreases. Goldsmith (1971) 
and Dasgupta (1986) have tried to explain this. It may be said to 
be due to the deformation of PEG molecules at pore entrance due to 
high presure.This deformation is also aggravated at higher 
concentrations . 

For BL the observed rejection increases with pressure 
differential initially and at high pressures it forms a 
plateau.This has also been confirmed earlier (Rajnish,et 
a.1 . , 1985;Poddar ,«•( a.1 . , 1989 ) . 

Though the actual reason for this relationship is not yet 
developed, one possible reason could be that BL forms a dynami c 
membrane (Goldsmith, 1971 ) .Since BL has a wide range of molecular 
weight , there is a possibility of a dynamic layer being formed by 
the bigger molecules .The plateau at higher pressure may be 
explained to be because of the characteristic of the layer 
formed. This dynamic layer may possibly be a gel layer with 
essentially constant membrane surface con centrtion.Hence, initial ly 
permeate concentration increases giving lower rejection and later 
stages shows a plateau. 
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Fig. 5-4 Effect of pressure differential on observed 
rejection at co = 47*12 rad/sec 






^ ^ - 2 Effect o£ Stirrer Speed on Observed Rejection : 

The rejection increases for both PEG and BL with increasing 
stirrer speeds. However ,this rejection attains constancy at high 
stirrer speeds. This may be explained as follows. 

As stirrer s peed increases , the feed solution becomes more 
uniform. This results in a lower membrane surface 
concentration, hence lower permeate concentration gives increasing 
rejection. A stage is reached where the solution is totally uniform 
an * no further effect of stirrer speed is noticed or observed 
rej ect ion . 

5.4.3 Ef f ect of Bulk Concentration on Observed Rejection s 

As evident from Fig. 5. 4 concentration always tends to 
increase the observed rejction fo BL.This may be due to the higher 
bulk concentration in Eqn.3.20,and also may be due to a marginal 
increase in permeate concentrations. Further, higher the bulk 
con cent rat ions, mo re are the possibilities of thicker layer 
deposited on membrane surface giving less amount of permeate 
solutes;hence teSs observed rejections. 

5.5 Fitting of Experimental Data : 

Once the datet of flux versus time was obtained, the next 
step would be to suitably represent the data by an equation with 
adjustable parameters .The model should be non-linear since the 
nature of flux versus time curve is highly non-linear. Due to the 
marked difference of curves (Figs. 5.1 and 5.2) at different set of 
experimental conditions, this difference has to be taken into 
account in the model equation. 

In favour of the above requirments the following model 
equation was chosen to fit the data! 



U 

J U f (t) * « + ft r X (5.1) 

In this equation a., ft and y are adjustable parameters , where y 
lies between O and 1 .The following are the advantages of choosing 
the model: 


(i) 

At 

time 

t = Of 

J uf (°) - « 

(ii) 

At 

time 

t )oo * 

J uf { ~) - « 


(since y lies between O and 1, the second term on the right of 
Eqn.5.1 is negligble with respect to the first term). 

( iii ) Although simple , the model was versatile and could fit 
any data within the required range of accuracy. 

(iv) Another most important advantage was that by measuring 
the zero flux,J uf (0) , the final flux J u ^,(oc ) and the UF flux J uf (t ) 
at any time (say t ^ ) one can calculate,}' .from Eqns.5. 1 ,5.2 and 
5. 3. This implies that the whole experiment could be directly 
modelled by just obtaining the data at three different 
points , provided the experimental data is accurate enough. 

(v) The model also satisfies the inherent nature of UF flux 
(Figs. 5. 1 ,5.2 and 5.3) i.e.,the drop in flux is ^ in an UF run^> 
falling with respect to time. 

Due to these advantages the model (Eqn.5.1) was fitted to 
each and every experiment and the parameters were estimated. 

To check the val i idi ty of the model , f or each run after 
estimating the parameters , an F-Test was conducted at a confidence 
level of 0.95. The computer program and the results are presented 
in appendix. 3 - 

The predictions of the model for all experiments is shown by 
a continuous bold line for BL and a continuous broken line for PEG 
in Figs. 5. 1,5.2 and 5. 3. The points shown in the same figures are 
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re&l experimental data. There exiata an excellent agreement between 
the experimental data and the model predictions - 
5.6 Analysis of Experimental Results ? 

5.6.1 Analysis with PEG-6000 : 

The deviation of the UF flux (Fig. 5. 3) from the solvent flux 
may be explained by either osmotic pressure model or gel layer 
model .The f ormer is more likely for lower molecular weight 

solutes ,whl ie the later for higher molecular solutes (Doshi , 1985) . 

PEG— 6000 is a lower molecular weight polymer. Also the 

molecular average weight of BL was found to be around 1500, since 
it was a degraded BL.So flux limitations through osmotic pressure 
model is most likely (Eqn.3.5). 

For a time dependent case Eqn.3.6 is a possible way to 

explain flux limitation. In Eqn .3.7, (F) was obtained from 

Eqn.5.1 at a constant flux stage with respect to time, which also 
corresponds to the final fluxes in Figs. 5 . 1 , 5.2 and 5 . 3. Further , it 
was observed that beyond 110 minutes the final flux ,J uf (F),for 

all the experiments attained the 1 imi t ing flux valuejhence this 

time was taken as the limiting time.R p (F) was then calculated from 
Eqn .3.7. A detailed procedure for calculating R p (F> is presented in 
the appendix (A— 4 ). The results are shown in Tables A— 4 . 1 . 

5-6.2 Analysis with Black Liq uor ? 

R (F) for BL was estimated in an identical way as for 
P 

PEG, except for the difference in calculating mass transfer 
coefficient k .Since k cannot be obtained using Eqn. 3. 10 as the 
properties of BL are unknown, velocity variation technique 
(Sec. 5. 3) was employed to estimate k .Equation 3.24 was plotted 
using experimental data for BL and the slope was obtained from 
Fig. 5. 5 -The results are shown in Table A-4.2. 
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5.7 Development of Correlation between R ^(F) and the System 
variables ' 

In UF systems it becomes necassary to have some idea of how 
the polarised layer resistance varies with the influencing 
parameters . In this work an attempt has been made to correlate 
R (F) with the operating variables. 

p 

R is a function of membrane surface concentration, flow 
P 

pattern (characterised by stirrer Reynold's No. ) , distance between 
the stirrer and the membrane surface, system geometry 
(baffles, channels, etc. , ) ,Froude No., applied pressure 

differential, etc. 

Since the geometry of the system is fixed for all the 

experiments, its effect can be neglected. Hence the effective 

variables are AP ,C and N .The stirrer Reynold’s No. is 

appl m Ke 

2 

defined by D Np/m - 

N c and AP , have been employed in the form of 
Re m appl 

dimensionless variables of the order of unity in order to obtain a 
simplified correlat ion.R p was made dimensionless by dividing by 

R +R .So the final form of the relationship becomes, 

sw a 


R (F) 
P 


R 


m 


+ R 


a 



AP 


appl 


An 
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The experimental values of 
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and N_ seperately . (Figs. 5. 6, 5. 7 and 5. 8). The 
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linearity obtained on log-log plots (Figs. 5. 6 and 5.7) implied a 
relationship of the forms 


R (F) 

P 

- = a r 

C 

m 

b 

) f 

AP 

appl 

C 

d 

R + R 
m a 

3 l 

p s 

J l 

An 

) 
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The values of a,b,c and d were estimated by least-square 
curve fit in a log-log form. After solving Eqn 5.5 the results 
obtained were , as follows: 

PEG-6000: 



The positive values for 'a'/b’ and*c* and negative value 
for *d* clearly indicates that the polarised layer resistance 
increases with C n /p % and AP^^^j/ATl , while it decreases with N^.The 
experimental and predicted values were plotted and are shown in 
Fig. 5. 9. The figure shows a good agreement between the correlation 
and the results. 

The deviations of BL may be due to pore clogging which 
increases the resistance, while it may also be due to excess 
passage of lower molecular weights of BL thus decreasing the 

resistance. 

The correlations (Eqns.5.6 and 5.7) may serve as an useful 
tool for estimating the extent of polarisation as a function of 





operating varlabl ta . Further , this may also predict the limiting 
flux phenomena to design an UF plant .However , the results obtained 
are for batch cell operated in a continuos mode. Hence further work 
is needed to generate experimental data out of tubular or spiral 
wound modules. 

5.8 Overall Analysis of Results ? 

A deep study of the Figs. 5. 6, 5. 7 and 5,8 coupled together 
with Figs .5. 1 ,5. 2, 5. 3 and 5.4 has to be done. 

The linearly increasing nature of R (F)/(R +R ) with respect 

p ' m a 

to C m /P s and AP^^/ATI and linearly decreasing nature with N^ e 
(Figs. 5. 6 and 5.7) suggests some conclusions. 

An increased drop in rejection (Fig. 5. 4) along with 
increasing R^(F)/ ( R m +R a > for PEG at high pressures suggests a 
boundary layer made of deformed PEG molecules but still resistant 
to the flow. 

For BL a constancy in rejection at high pressures (Fig. 5. 4) 
with a late limitation (Fig. 5.1) and increasing R^(F) (Fig. 5. 7) 
suggests continuous formation of polarised layer. This may be 
elaborated that as the pressure increases convective flux continue 
to increase which continue to form polarised layer on the membrane 
surface. This phenomena exihibits continuous fall of flux as the 
pressure increases. 

(iii)Higher stirrer speed forms an early limitation 
(Fig. 5. 2) , higher rejection (Fig. 5. 4) and lesser R p <F) (Figs. 5. 7 
and 5.8) for BL.This suggests a lower back diffusive flux. So in 
this case osmotic pressure may not be controlling. 



CHAPTER-6 


CONCLUSIONS AND RECOMMENDATIONS 

Ultrafiltration of PEG-6000 and Black liquor solutions using 
an asymmetric membrane and analysis of the experiment^ results 
as was done in this work leads to the following conclusions: 

(i) Black liquor has a significant tendency to form a 
polarised layer while PEG does not have so. 

( ii ) Continuous formation of polarised layer on membrane 
surface with increasing pressure results in a continuous fall of 
permeate flux for BL while insufficient formation of polarised 
layer for PEG does not exhibit significant fall in permeate flux 
as a function of time. 

(iii) Higher convective flux due to higher Reynold’s 
No., tends to precipitate BL on membrane surface earlier showing 
limitations in permeate flux compared to lower Reynold’s No. 

Civ) PEG solutions behaviour may be concluded by the control 
of osmotic pressure model jwhile BL “ - ^ H- l3 r 

• r 1 is. not obeying totally osmotic pressure 

model .Hence, polarised layer deposition phenomena t ruely represents 
the limiting flux characteristics of BL. 

The following recommendations are made for future work in 
this area: 

( i ) Intermediate values of R (t) should be found in order 

P 

to detect the exact nature of variation of the polarised layer 


resistance. 





(ii) Th« mb.bb transfer coefficient has to be deternined at 
the membrane surface concentr at i on • To achieve this an i te rat i ve 
way could be adopted to find the dif f usi vi t ies .This also needs the 
physical properties of PEG and BL to be estimated beforehand at 
such high concentrations. 

( i i i ) The thoer ies of adsorpt ion or interact ion of solute to 

membrane should be investigated to quantify adsorption resistance 

(R ) • 
a 
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APPENDIX A 


A— 1 Determination of Viscosity— Concentration Correlation for 


Black-Liauor and PEG-6000. 



TABLE A— 1 

■ 1 Viscosity- 

-Concentration 

Data for PEG-6000 

So lutes PEG-6000 



Temperature 30‘C 




C t 
s s 

Concentration 

Time of Fall 

Density 

P s^ 0. 999567t . 

w 

X 

t(sec) 

gm/cc 

cP 



* 

* 

0 (water) 

360.6 

0.995678 

0.8500000 

1.0 

431.67 

0.99809963 

1.0200000 

2.5 

561.40 

1.0020412 

1.3317794 

5.0 

801.20 

1.0068429 

1.9097190 

10 

1506.4 

1.0136519 

3.6149594 


* Handbook for Chemical Engineers, Perry, XH. , VI Edn, 


McGraw-Hi 11, NewYork. 





TABLE A-l . 2 Vlacoalty-Coneentrat ion Data, for BL 


SolutesBL Temperatures30"C 


Concent rat ion 

X 

Time of Fall 

t (sec) 

Density 

gm/cc 

C t 

s s 

Hs 0.995678t 
cP W 



*§■ 

* 

0(water ) 

360.60000 

0.9956780 

0.850000 

in 

a 

o 

371.76800 

0.9977921 

0.878188 

1.0 

381 . 14480 

0.999414 

0.902275 

2.0 

396.46000 

1.005848 

0.944088 

4.495 

447.0660 

1.020878 

1.080488 


The following correlations have been obtained for the 
viscosity-concentration data by polynomial fit of the 
data(Poddar,T.K. , 1986) . 

2 -33 -44 

^PEG “ 0,85 + °- 1446 C + 0-02734 C -4.276x10 C +2.84x10 C 

( A-1 ) 

—2 2 

u , = 0.85 + 0.06147 C - 1.117x10 C + 1.98x10~ 3 C 3 (A-2) 

Bl_ 

where C is the concentrat ion of the solution expressed in percent 
(X) gm of solute in lOOcc of solution. 

* Handbook for Chemical Engineers , Perry .J?.H, VI Edn, 


McGraw Hill ,NewYork 



A- 2 Determination of Denaity of Solution for Black Liquor and 
PEG-6000 s 


TABLE A-g Density-Concentration Data for PEG and BL 


Solution 

Concentration 

{•/.) 

Specific Gravity 

Density 

gm/cc 


1.0 

1 . 00E4E88 

0.9980963 


E.5 

1.0063908 

1.00E041E 

PEG-6000 

5.0 

1 .01 1E134 

1 . 00684E9 


10.0 

1.0180519 

1.0136519 


0.5 

1.00E1EE3 

0.99779E1 


1.0 

1 . 004E8E0 

0.9999414 

BL 

E.O 

1 . 010E143 

1 . 005848E 


4.495 

1 . 010E143 

1.0EE878 


A— 3 Determination of Concentration of PEG-6000 by Ref ractometer ? 

Refractive Index of PEG-6000 at various concentrations were 
measured at 30*C.The Fig.A-1 shows the calibration curve of 
Refractive Index versus concentration (gm per lOOcc) of 
PEG-6000. Table A-3 shows the Refractive index-concentration data. 
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Fig. A-1 Calibration curve for PEG-6000 concentration 
measurement 



TABLE A-3 Refractlvelndex-Concentratlon Data for PEG 


Temperature: 30 *C 


Concentration 

X 

R.I. 

Concentration 

R.I 

O 

» 

O 

1.3317 

2.0 

1.3345 

0.2 

1.3320 

3.0 

1.3360 

0.5 

1.3325 

4.0 

1.3376 

0.8 

1.3329 

5.0 

1.3381 

1.0 

1.3331 




A— 4 Calculation of Flux, Observed Rejection and R (F> : 

The following steps were carried out in a sequential order 
to calculate various quantities to estimate R p (F) . 

(i) The R value was obtained from Eqn.3.1 using the pure 

m 

solvent flux for the fresh membrane. 

(ii) The value for R +R of after every experiment was 

m a 

calculated since the pure solvent flux was already obtained from 
experiment (Eqn. 3. 3) . 

(iii) For a paritcular feed concentration (C^) ,mass transfer 

coefficient (k) was calculated from Eqns.3.10 and 3.11 for PEG. For 
Black Liquor the velocity variation technique (Sec. 3. 5 and 

Fig. 5. 5) was employed to calculate the mass transf er coefficient. 
For these calculations the physical properties were obtained from 
Tables A— 1 and A— 2. 

(iv) The flux of the ultrafiltration experiments was fitted 


to model equation (Eqn .5.1). 
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(v) The final flux J u ^(F) vaa calculated from Eqn.5.1 at a 
time T(F) (110 min) when there was no change in UF flux 
(Sec. 5. 6.1). 

(vi) The membrane surface concentration (C ) was then 
calculated from Eqn.3.9 using the final flux J u ^(F) from step v. 

(vii) The osmotic pressures for the permeate (n^J^and for the 

membrane surface (n ) were calculated from Eqns.3.15 through 3.19 

m 

using the experimental value of C p and the calculated value of C m - 
Subsequently An !«, - n p ) was calculated. 

(viii) R (F) was then calculated , since all other quantities 
P 

are known in Eqn 3.7. 

(ix) The above steps were repeated for all runs. 



TABLE A- 4.1 


PEG-6000 


EXP AP 

C b 

to 

C 

P 


P 


REFNO N/m 2 

kg/m 2 

rad/s 

kg/m 2 

kg/ms 

xIO 3 

kg/m 3 
xIO 3 

m 2 s 

10 6 

00 

55163.3 

16.667 

20.944 

1.333 

1 . 1527 

0.9985 

1 . 1528 

01 

55163.3 

16.667 

47.124 

2.750 

1 . 1527 

0.9985 

1.1528 

02 

55163.3 

9.167 

47.124 

0.2216 

1.0038 

0.99809 

1 . 0057 

03 

55163.3 

44.670 

47.124 

3.5336 

1.7823 

1.0058 

1.7720 

04 

55163.3 

16.667 

34.034 

3.0000 

1.1526 

0.9999 

1 . 1528 

05 

55163.3 

16.000 

5.502 

1.3330 

1 . 1388 

0.9967 

1.1426 

06 

689516.6 

16.667 

47.124 

2.000 

1 . 15266 

0.9985 

1.1528 

07 

275806.6 

16.667 

47.124 

2.000 

1 . 15266 

0.9985 

1.1528 

08 

55163.3 

44.67 

34.034 

7.857 

1.7823 

1.0058 

1.7720 

09 

55163.3 

43.83 

55.502 

3.803 

1.7625 

1.0057 

1.7526 

10 

689516.6 

41.33 

47.124 

8.333 

1.704 

1.0051 

1.6951 

11 

275806.0 

44.00 

47.124 

2.0 

1.7665 

1.0056 

1.7565 


continued 




Table A-4.1 (continued) 



k 

xIO 5 

m/s 

J ,<F) 
uf 5 

xIO 

m/s 

C 

m 

kg/m 3 

X 

r 

1m 

r 

2m 

2 

xIO 

r 

Ip 

r 

Hp 

xIO 3 

1.009 

1.4742 

67.410 

331.72 

0.940 

5.992 

0 , 9988 

1.184 

1.931 

1.8899 

39.791 

331.72 

0.9646 

3.537 

0.9975 

2.44 

2.058 

2.435 

29.364 

332.29 

0.965 

3.499 

0.9980 

1.967 

1.579 

1.068 

84.413 

329.75 

0.9249 

7.503 

0.9968 

3.140 

1.488 

1.326 

36.322 

331.72 

0.9677 

3.229 

0.9973 

2.666 

2.210 

2.045 

38.331 

331.72 

0.9659 

3.401 

0.9988 

1.184 

1.931 

2.072 

44.899 

331.72 

0.9600 

3.990 

0.9982 

1.777 

1.931 

0.935 

25.813 

331.72 

0.9770 

2.229 

0.9782 

1 .777 

1.217 

0.751 

76.03 

329.29 

0.9324 

6.758 

0.993 

6.894 

1.810 

1.116 

77.908 

329.8 

0.9307 

6.925 

0.9966 

3.380 

1.615 

1.169 

76.362 

329.95 

0.9321 

6.787 

0.9925 

7.407 

1.586 

5.569 

61.665 

329.79 

0.9451 

5.481 

0.9982 

1.777 


. . . continued 
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Table A-4.1 (continued) 


n 

m 

N/m 2 

n 

P 

N/m 2 

An 

N/m 2 

J'kIO 5 

s 

m/s 

R +R xIO 
m a 

ki 2/ 2. 

Ns /m kg 

13 _ . -13 

R (F)x10 

P 

Ns 2 /m 2 kg 

60758.5 

510.49 

60248 . 1 

3.5178 

1.8447 

1 . 1321 

25741.6 

1152.18 

24589.4 

3.5415 

1.8324 

0.5744 

25373.6 

86.018 

25287.6 

3.7835 

1.7152 

0.3246 

91818.8 

1400.7 

90418.1 

3.7502 

1.7304 

0.7072 

22427 . 1 

1225.6 

21201.5 

3.6427 

1.7814 

1.7393 

32124.4 

1898.2 

30225.8 

3.6726 

1.7670 

0.4571 

30895.1 

810.3 

30084.7 

3.7622 

1.7249 

0.9941 

13882.4 

810.3 

13072.1 

3.5234 

1.8418 

0.7116 

75491.2 

3308 . 8 

72182.3 

3.7039 

1.7520 

1.7426 

78911.1 

1517.0 

77394.0 

3.7304 

1.7396 

0.6583 

76038.1 

3539.6 

72498.5 

3.7529 

1.7292 

1.3583 

52051.8 

786.8 

51264.9 

3.8159 

1.700 

0.5814 


. . . continued 
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Table A-4.1 (continued) 


R (F) 

P 

C m 

AP , 

appl 

N Re 

R +R 
m a 


An 

0.643433 

.006742 

9.15569 

6387 11 

0.31141 

0.03979 

22.4329 

14371.0 

0.17718 

0.02936 

21.8135 

16500.51 

0.4 1232 

0.08392 

6. 10069 

9349.42 

1.0051 

0.03632 

26.0176 

6387.05 

0.25658 

0.03834 

18.2497 

16965.9 

0.56253 

0.0449 

22.9196 

14371.0 

0.041257 

0.02581 

21.09887 

14371.0 

0.94615 

0.07559 

7.6419 

6752.36 

0.37572 

0.07747 

7.12733 

11028.4 

0.7808 

0.07596 

9.5107 

9773.63 

0.33624 

0.06131 

5.38002 

9431.86 


TABLE A-4.2 


BLACK LIQUOR 


EXP AP 

C b 

63 

C 

P 


P 


REFNO N/m 2 

kg/m 2 

rad/ s 

kg/m 2 

kg/ms 

xIO 3 

kg/m 2 

xIO -3 

2 

m s 

10 6 

20 

55163.3 

4.083 

47.12 

2.845 

0.8733 

0.9977 

0.8753 

21 

55163.3 

8.194 

47.12 

4.363 

0.8939 

0.9991 

0.8947 

22 

275806.6 

17.249 

47.12 

9.947 

0.9329 

1.004 

0.9290 

23 

55163.3 

16.317 

47.12 

8.888 

0.9291 

1.003 

0.9257 

24 

689516.6 

27.093 

47.12 

9.982 

0.9739 

1.010 

0.9641 

25 

55163.3 

17.228 

34.03 

9.1012 

0.932 

1.004 

0.9289 

26 

55163.3 

18.082 

55.50 

9.327 

0.9364 

1.004 

0.9319 

27 

55163.3 

44.806 

47.12 

25.016 

1.0792 

1.0207 

1.0572 

28 

55163.3 

43.393 

34.03 

23.32 

1.0681 

1.0199 

1.0473 


continued 





Table A-4.2 (continued) 


J -(F) 
uf«- 

xIO 

m/s 

R o 

Y. 

k 

xIO 5 

m/s 

C 

m 

kg/m 3 

X 

r 

1m 

r 

2m 

2 

xIO 

r 

Ip 

r 

2p 

xIO^ 

2.988 

30.2 

4.768 

5.161 

392.7 

0.995 

0.4588 

0.9974 

2.528 

2.45 

46.57 

4.768 

10.76 

392.2 

0.990 

0.957 

0.9961 

3.878 

1.015 

42.39 

4.768 

19.25 

390.27 

0.983 

1.711 

0.991 

8.841 

1.809 

45.52 

4.768 

19.74 

390.48 

0.982 

1.756 

0.992 

7.900 

2.334 

63.15 

4.768 

37.90 

387.99 

0.966 

3.369 

0.991 

8.873 

1.281 

47.17 

3.675 

20.61 

390.27 

0.981 

1.832 

0.991 

8.09 

2.107 

48.41 

5.43 

22.94 

390.08 

0.979 

2.039 

0.992 

8.291 

1.041 

55.83 

0.644 

49.63 

383.94 

0.955 

4.412 

0.9778 

22.23 

0.871 

46.26 

0.497 

47.41 

384.25 

0.957 

4.214 

0.979 

20.73 


continued 
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Table A-4.2 (continued) 


n 

m 

N/m 2 

n 

P 

N/m 2 

An 

N/m 2 

J'xlO 5 

s 

m/ s 

— 

R +R xIO 
m a 

Ns 2 /m 2 kg 

R (F>x10 

P 

k , 2. 2. 

Ns /m kg 

17P1.93 

944.64 

837.283 

4.0285 

1.6108 

0.41383 

4076.64 

1495.45 

2581.19 

3.9097 

1.6598 

0.8957 

8395.41 

3735.67 

4659.75 

3.7829 

1.7154 

1.2017 

8657.82 

3280.75 

5387.06 

3.6173 

1.794 

1.53418 

21830.8 

3779.58 

18051.3 

3.3149 

1.9576 

1 . 1589 

9182.09 

3160.55 

6021.54 

3.509 

1.8493 

2.6059 

10593.2 

3479.09 

7114.14 

3.429 

1.8925 

0.9098 

33698.2 

12041.5 

21656.7 

3.362 

1.893 

2.8198 

31246.1 

10937.9 

20308.2 

3.3672 

1.9272 

3.8173 


continued 
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Table A-4.2 (continued) 


R (F> 

P 

C * 

AP , 

appl 

N Re 

R +R 

m m 


ah 

0.2439 

0.00516 

658.81 

18927.68 

0.5560 

0.01077 

213.70 

18517.27 

0.7239 

0.01917 

59.189 

17833.05 

0.8943 

0.01967 

102.39 

17895.96 

0.6459 

0.03752 

38.197 

17183.21 

1.3311 

0.0205 

114.51 

12880.44 

0.4919 

0.02284 

91.606 

20741.33 

1.4899 

0.04863 

25.470 

15669.67 

2.0165 

0.04648 

27.162 

11424.94 
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************************************************************** 
REFERENCES: TMSL-MANUAL. , HANDBOOK OF CHEMICAL ENGINEERS } 

PERRY (EUR ) , VI EDN.TI-CHAPTER(STATI5TTC5) .MCGRAW HILL. 
NOMENCLATURE : ESS S ERROR SUM OF SQUARES ,RSS : REGRESSION 
SUM OF SOUARES.MSRSSI MEAN SQUARE OF REGRESSION SUM OF 
SQUARES.MSESSjMEAN SQUARE Of ERROR SUM OF SQUARES 
SEE: STANDARD FRROP ESTIMATE, RSMTT7.:TMSL ROUTINE 
ESAMPLlF ESTTMATF OF DATA 

*******************41****************************************** 
DIMENSION XY(9,J1 »TVf i) ,PAR£3) ,STAT(6) , YPREDf9) . 

TTERf 2) ,VCV(6,6) 

EXTERNAL RSMSSE . UERTST . ZXFTB 
REAL MSRSS.MSESS 

OPEN . ( nNTT»33.DEVlCE='nsK' , FT LFs 'EY PO 1 .100') 

OPEN r UMIT=6, DEV ICES 'DSK') 

OPEN ( UNTT=32. DEVICE= 'DSK ',FILFr 'EXPO 1. OUT') 

T V f 1 ) =9 

READf33,*)fXY(I,l),XY(T,2),I=l,iV(l)) 

,******************$******* ********* **************************** 
IXYsIVfl) 

T V ( 2 ) =2 
T V ( 3 ) =0 
PAR ( 1 ) =0^0 
PAR (2) = 1„0 
PAR ( 3 ) =0 . 0001 
PAUSE 1 . 

CALL RSMTT7.(XY,IXY»IV, PAR. STAT.VCV. ITER. TER) - 

**************** ANOVA ****¥**»*♦***■*************************** 
ESSsSTATf 5 ) 

DF=IVC1 )~2 

MSESS=ESS/DF 

RSS=0 

DO !“> Lsl , T V f 1 ) 

YPRED(L)sSTAT( 1 )+STAT( 7)* (STATC 3) ♦♦XY(L, 1 ) ) 

RSS = RSS+f YPREUf L)-STATf4) )**?. 

CONTINUE 
MSRSS=RSS/1 .0 
ESAMPL=MSRSS/MSESS 
SEEsSURTrESS/fTVf 1 )-3 ) )' 

******************************************* 3 ,.^**^*******^******* 

sssssssssssssssssssssssssssssssssssssssssssssssssssssssssssssss 

WRITE ( 32 . 100 ) 

WRITE ( 32.98) STA^Cl) 

WRTTE(32.97)STAT(2) 

WRTTE(32.96)STATf 3) 

WH,TTE(32.95)STATf4) 

WRITE ( 32 . 94) STAT ( 5 ) 

WRTTE£32.93)STAT(6) 

WRITE (32. 91) 

WRITE (32.98) 

WRTTE (32.92) (XY(M, 1 
WR.TTE( 32.89) 

WRTTE(32.90) ((VCV(M.L) ,M=1,3) ,L=1 
WRITE C5.*)IER 
WRTTE(5,*) ITER 
TYPE* , ES AMPL 

F0RMATC//15X. 'RESULTS OF. REGRESSION FoR EXP01 '/) 

F0RMAT(//5X 'ALpHA=' .E14.8) 

E0RMATC5X, 'BET As ' . F 1 4 . 8 ) 

F0RMAT(5X, 'GAMMa=',E14.8) 

E0RMAT(5X, 'AVG OF DEP VAR Ys'.ElA.R) 

FORMAT (5X, 'ERR SUM SOUARESs ' , El 4. 8 ) . 

F0RMAT(5X. 'AVG OF GAMMA POWER Xs'.El4„8). 

F0RMAT(5X.E14 8 # 5X f El 4.8.6X.E14.8.5X-.E14.8) • 

FORM AT (/'*******♦*** ************************«**********') 
EORMATf/lOX'SAMpT.E VARIANCE-CqVaR A lNCE MATRIX'/) 
F0RMAT(10X,E14.8) .. , 

E0RMATC5X. 'INDEP VArTABLF' , 5X . 'DEPEND VARIABT.E'. 

5X 'GAMA . POWER X ' , 8 X . ' YPR ED ' / ) 

WRITE (32,88) 

WRTTE£32.87)ESS.MSESS 
WRITE £32.86)RSS ,'MSRSS , FSAMPL 
WRITE £32. 84) SEE 
WRTTE ( 32 . 83 ) 

EORMATf/RX.'ERR SS= ' . El 4 . 8 , 5X . ' ME AN Sqs'F!4.R) 

FORMAT (/'*************** ANOVA AND F-EXpOl ******************** 
FORMAT(/?X.'REG SS='.E14.8.5X.'MEaN SQs 'E l 4 . ft . 5X . 

'F SAMPLEs',Fl0.5) 

FORM AT (/7X. 'STANDARD ERROR, ESTIMATES '. El 4 . R ) . 

FORM AT (/'*************************************************') 

STOP 

END 


. XY ( M # 2 ) ,XY (M,3) , yPRED ( M) , Msl , IV ( 1 ) 3 
3) 



RESULTS OF REGRESSION FOR EXPOO RESULTS OF REGRESSION FOR EXP01 
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